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Abstract 

Angle-resolved photoemission spectroscopy (ARPES) is a powerful experimental technique for 
directly probing electron dynamics in solids. The energy vs. momentum dispersion relations and the 
associated spectral broadenings measured by ARPES provide a wealth of information on quantum 
many-body interaction effects. In particular, ARPES allows studies of the Coulomb interaction 
among electrons (electron-electron interactions) and the interaction between electrons and lattice 
vibrations (electron-phonon interactions). Here, we report ah initio simulations of the ARPES 
spectra of graphene including both electron-electron and electron-phonon interactions on the same 
footing. Our calculations reproduce some of the key experimental observations related to many- 
body effects, including the indication of a mismatch between the upper and lower halves of the 
Dirac cone. 
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In ARPES experiments a sample is illuminated by monochromatic photons, which can 
extract electrons from the sample if the photon energy exceeds the work function. Analysis 
of the kinetic energy and angular distribution of the emitted electrons yields the binding 
energy of the electron in the material and its crystal momentum parallel to the surface^. 
The measured intensity /(k, cj), where k and u are the momentum and the binding energy 
of electrons (usually referenced to the Fermi energy), can be written as^ 

/(k, w, hv) = /o(k, u, e^, hu)f{uj)A(k, u) , (1) 

where the function /o(k, w, Cj^, /iz/) takes into account the absorption cross section of the 
incident photon of energy hu and polarization e^. The function f{uj) is the Fermi-Dirac 
distribution, and v4(k, u) is the electron spectral function^. In interpreting ARPES mea- 
surements in a narrow energy range, it is appropriate to assume that the absorption cross 
section of the photon is constant. Under these conditions, a measurement of the ARPES 
spectrum provides direct access to the electronic spectral function v4(k, u). Within quantum 
many-body theory the electronic spectral function can be expressed as 



2 -ImS(k,a;) 
vr [w - £k - ReS(k, uj)f + [ImS(k, uj)f 



A{K^) = Z1 ^ rx ^12 > (2) 



where the ^k's are the the single-particle energy eigenvalues of a reference mean-field system, 
and the self-energy E(k, u) accounts for the many-body interactions going beyond the mean- 
field picture^. For simplicity of discussion, the band indices are dropped from Eq. ([2]) (see 
Methods). 

In our investigation, we use density-functional Kohn-Sham eigenstates to describe the 
mean- field or non- interacting electrons. The electron self-energy arising from the e-e interac- 
tion T,^~'^(k,Lo) is evaluated within the GqIVq approximation^ (Fig. la). This corresponds to 
retaining the first diagram in the Feynman-Dyson perturbation expansion of the self-energy 
operator but in terms of the screened Coulomb interaction Wq (Supplementary Discussion 3 
for a comparison with previous first-principles calculations^!^). In this work, Gq is con- 
structed from the Kohn-Sham eigenvalues and eigenfunctions of density-functional theory, 
and Wo is the bare Coulomb interaction screened by the full frequency-dependent dielec- 
tric matrix e(r,r',a;) calculated within the random phase approximation. The self-energy 
S'^^P^(k, cj) arising from the e-ph interaction is similarly evaluated within the Migdal ap- 
proximation^ (Fig- lb). With these choices, the e-e interaction and the e-ph interaction are 



described consistently within the same level of approximation^ (see Methods). The total self 
energy is then obtained as 



Graphene, a single layer of carbon atoms in a honeycomb structure, has recently become 
an active research area in physics, chemistry, and nanoscience not only because of its peculiar 
low-energy massless Dirac fermion band structure^"-, but also because it holds promise 
for novel electronics and spintronics applications^^. In particular, the epitaxial growth of 
graphene on silicon carbide (SiC) has emerged as one of the promising routes towards large- 
scale production of graphene^^i^. 

The interpretation of the measured ARPES spectra of epitaxial graphene grown on silicon- 
rich surface of SiC has been controversial. The spectral features observed in early ARPES 
measurements^ were interpreted qualitatively in terms of e-e and e-ph interactions. On 
the other hand, experiments and analyses performed by a different group suggested that 
the low-energy ARPES spectrum is dominated by a quasiparticle energy gap of 0.2-0.3 eV 
at the Dirac point^^. According to Ref.-i^, this band gap likely arises from the coupling of 
the graphene layer with the reconstructed surface of the SiC substrate. Despite a number 
of subsequent studies to resolve this controversy^i^^, the detailed nature of the low-energy 
quasiparticle dynamics in epitaxial graphene remains an open question. 

The electronic structure and the photoemission spectra of graphene have also been ex- 
plored in a number of theoretical investigations using density-functional theory approaches 
within a non-interacting single-particle picturei^^-ii. The effects of many-body interactions 
have also been investigated in model calculations^i^^. However, first-principles calculations 
of the full k- and w-dependent ARPES spectral function - which contains both the quasi- 
particle dispersions and their lifetimes - including the e-e and the e-ph interactions have not 
been reported. The lack of first-principles many-body investigations can partly be ascribed 
to the extremely demanding computational efforts required to evaluate the real part of the 
electron self-energy, both for the e-e and for the e-ph contributions. 

In dealing with effects of e-e interactions, first-principles calculations have advantages 
over model calculations based on the two-dimensional massless Dirac equation^i^^. First, 
the scattering rate of even the low-energy charge carriers, whose non-interacting dispersion 
relation can be well approximated by the two-dimensional massless Dirac equation, shows 




(3) 
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strong wavevector- anisotropic behaviors^°. This is because the carrier scattering rate in 
graphene depends sensitively on the sign of the band curvature^^. Therefore, the significant 
wavevector anisotropy in the electron scattering rate is not captured by the model calcu- 
lations. Second, the model calculations require a cutoff of the high-energy states and an 
adjustable parameter mimicking effects of internal screening arising from the high-energy 
states (including the vr states higher in energy than the cutoff and the a states), in addition to 
the external screening due to the environment. The first-principles approach employing the 
full bandstructure accounts for these processes explicitly and requires neither a high-energy 
cutoff nor an empirical parameter to describe the internal screening. 

To determine the quasiparticle energy vs. momentum dispersion relations from the cal- 
culated ARPES intensity maps for graphene (Figs. 2a, 2d, and 2g), we follow the standard 
procedure adopted in analyzing ARPES experiments^. First, the energy distribution curves 
(EDCs) are obtained by performing cuts of the intensity maps at fixed photoelectron momen- 
tum (vertical cuts in Figs. 2a, 2d, and 2g). Subsequently, the quasiparticle band structures 
are generated by connecting the locations of the maxima in the EDCs for each photoelec- 
tron momentum (Figs. 2b, 2e, and 2h). This procedure ensures that the calculated and the 
measured dispersion relations are obtained from the corresponding intensity maps using the 
same procedure. 

While the dispersion relations extracted from our spectral functions for an isolated 
graphene layer are linear at large binding energy, we observe a sizable kink for n-doped 
graphene near the Dirac point (k = in Fig. 2) at an energy below the Dirac point energy. 
Such kinks result in a mismatch between the linear extrapolations of the lower and the 
upper portions of the Dirac cone (Fig. 2). A similar phenomenon has been observed in the 
measured ARPES spectrap^ii^. In order to quantify the size of the kinks and the associated 
energy mismatch, we have taken the energy difference between the two linear asymptotes 
of the upper and the lower bands (Akink in Fig. 2f). This calculated energy offset Akink is 
predominantly a result of many-body effects, which comes from the GVF self-energy, and is 
found to increase with the doping level (see Fig. 3), consistent with the experimental trenci^. 
We have checked that, if many-electron effects are not considered, the energy mismatch is 
several times smaller, with value < 20 meV for the most heavily doped case considered here 
and even smaller for other cases. We note that, for the path along which the ARPES spectra 
are calculated (inset of Fig. 2c), the nonlinearity of the bare graphene band is smallest. 
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Although the calculated Akink shows qualitative agreement with experiment, after taking 
into account the screening of the SiC substrate (Refs.— see Methods), our calculated 
values underestimate the experimentally observed offsets in Ref.— consistently by 60-90 meV; 
this discrepancy, together with the comparison between theory and experiment of electron 
linewidths discussed in Ref.~ and later in this paper, suggests the possible opening up 
of a band gap at the Dirac point energy due to the interaction between graphene and 
a reconstructed silicon-rich surface of SiCi^»i^""— . Our study shows the need of further 
first-principles studies, considering both the atomistic structure of the graphene-substrate 
interface and many-body effects. 

The quasiparticle velocity can be extracted from the simulated spectral functions. We 
find that e-e interactions greatly enhance the band velocity by over 30 % compared to the 
density-functional theory value in pristine graphene, but dielectric screening from the SiC 
substrate (by weakening the e-e interaction) reduces the quasiparticle velocity by as much as 
~10 % (Supplementary Discussion 1 and Fig. [5]). Moreover, the calculated velocity decreases 
with doping (Fig. [5]) in agreement with previous calculations^. 

Our calculations also reveal phonon-induced kinks near the Fermi energy at binding 
energies between 150 and 200 meV (e. g. , Fig. 4a) in good agreement with experimental 
photoemission maps^i"— . These signatures of the e-ph interaction in graphene have been 
analyzed thoroughly both experimentally and theoretically22i"— . 

A complementary and important piece of information provided by ARPES intensity maps 
is the linewidth of the electronic quasiparticle peaks. The linewidth r„k is related to the 
lifetime r„k of the electron in a given quasiparticle state through r^k = 2^/r„k, and plays 
an important role in transport phenomena. The electron linewidths are extracted from the 
measured photoemission spectra by considering momentum distribution curves (MDCs), 
which are constant-energy cuts of the intensity maps. The width of the MDC at a given 
binding energy can subsequently be obtained using a Lorentzian fit^. From our simulated 
ARPES spectra of graphene along the FKM direction in wavevector space (Fig. 4a), we 
obtained the MDCs (Fig. 4b) and their widths (Fig. 4c). Unlike previous methods^, this 
procedure is direct and is not confined to materials having a linear electronic dispersion. 

The widths extracted from the calculated MDCs in n-doped graphene after including 
the effect of the SiC substrate screening (see Methods) follow closely the experimental mea- 
surements on epitaxial graphene at large binding energies^, although they underestimate 
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the experimental data close to the Dirac point energy (Fig. 4c). The agreement between 
calculated and measured widths at large binding energy and the underestimation of the 
linewidths close to the Dirac point energy provide additional support to the proposed sce- 
nario of a band gap opening arising from the interaction with the SiC substrate^^ii^"— i^. 
That is, the opening of a band gap and the generation of midgap states near the Dirac 
point energy (not included in our calculations) would lead to increased linewidths in this 
energy regime^^. However, from the view point of theory, a conclusive statement can be 
made only after first-principles calculations considering both the atomistic structure of the 
graphene-substrate interface and the many-body effects, which is beyond the scope of this 
work. 

The widths extracted from the calculated full ARPES spectrum and those obtained by 
calculating the imaginary part of the on-shell electron self-energy Im S(k, Sk)^ are similar 
for n-doped graphene with the same charge density. (The charge density of n-doped graphene 
considered in Ref.— is different from that in Fig. HI) However, it is important to note that the 
current method of using the full ARPES spectra, although involving heavier computations, 
is more powerful because it can be in general applied to systems whose bare electronic energy 
dispersion is not linear. 

In this work, the effects of substrate optical phonons have not been considered. These 
substrate phonons contribute to the room-temperature transport properties of graphene^^. 
However, over a 1.5 eV range around the Fermi level in doped graphene relevant to ARPES 
experimental*^, the imaginary part of the electron self energy that arises from interactions 
with the substrate phonons is less than 1 meV^, about two orders of magnitude smaller 
than the self-energy effects intrinsic to graphene itself. 

The present work shows that first-principles simulation of ARPES spectra based on a 
quantum many-body theory approach treating e-e and e-ph interactions on the same footing 
holds great potential for the interpretation of complex ARPES spectra. In particular, a direct 
calculation of the quasiparticle spectral function is needed to obtain meaningful comparisons 
with experimental data (e. g. , the extraction of EDCs, MDCs, linewidths, and quasiparticle 
dispersion relations). 

Methods 
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1. Computational setup 

The Kohn-Sham eigenstates of graphene are obtained using density-functional theory calcu- 
lations within the local density approximation (LDA)^ in a supercell geometry2£. Electronic 
wavefunctions in a 72 x 72 x 1 k-grid are expanded in a plane-waves basis with a kinetic 
energy cutoff of 60 Ry. The core-valence interaction is treated by means of ab initio norm- 
conserving pseudopotentials^. Graphene layers are separated by 8.0 A and the Coulomb 
interaction is truncated to prevent spurious interaction between periodic replicas^. We 
have checked that increasing the interlayer distance to 16.0 A makes virtually no difference 
in the calculated self energy. Charge doping is modeled by an added electron density with 
a neutralizing background. 

Extending the procedure presented in Ref.-^, where only the imaginary part of the on-shell 
electron self-energy lmS^~'^{e-k) arising from electron-electron (e-e) interactions is calculated 
(^k being the LDA energy eigenvalue), we calculate the full frequency dependence of both 
the real and the imaginary parts of the dielectric matrix (within the random phase ap- 
proximation) and the self-energy operator Ti'^^'^{uj) within the GqWq approximation in the 
present work. Thus our theory includes the two scattering mechanisms arising from e-e 
interaction effects discussed in previous model calculations performed within the massless 
Dirac equation formalisro^^ii^, i. e. , electron-hole pair and plasmon excitations. Since our 
calculations are based on first-principles, they are parameter free for suspended graphene 
and can give information that depends on atomistic details, e. g. , the effects of trigonal 
warping. For convergence of the real part of T,'^~'^{uj), we have included conduction bands 
with kinetic energy up to 100 eV above the Fermi level. The frequency dependent dielec- 
tric matrix eG.G'(^)'^) is calculated within the random phase approximation using the LDA 
wavefunctions on a regular grid of u with spacing Au =0.125 eV^^, and the dielectric ma- 
trix at energies in between frequency grid points is obtained by a linear interpolation. In 
the calculation of the polarizability, for numerical convergence, an imaginary component 
of magnitude Au of 0.125 eV is introduced in the energy denominator. Convergence tests 
showed that the dimension of the dielectric matrix may be truncated at a kinetic energy 
cutoff of h'^G'^/2m =12 Ry. Additionally, we obtain the electron self-energy arising from 
electron-phonon (e-ph) interactions T,'^~^^{uj) following Ref.^^ for different levels of doping. 

In our approach, we describe the e-e and the e-ph interactions within the same level of 
approximation. Previous studies suggested the use of the quasiparticle dispersions renor- 
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malized by the e-e interaction to compute the e-ph interaction^^. This latter procedure (not 
adopted here) would correspond to including some of the higher order processes in ll^~^^{uj) 
whilst neglecting them in and would result in an unbalanced evaluation of e-e and 

e-ph effects according to different levels of approximation. In any event, even if some e-ph 
matrix elements calculated after renormalizing the bands through the GH^ approximation 
were 20% to 40% larger than those used in this work^^, the band velocity would also be 
enhanced by a similar factor (this work and Refs.-i^), and the two factors would cancel out 
approximately in the calculation of E'^^p'^(w) [cf. Eq. (2) of Ref.— ]. Therefore, we estimate 
that the effects of such alterations on our results are not significant. 

2. Angle-resolved photoemission spectra 

First, the trace of the spectral function with respect to band index n, i.e., yl(k, w) = 
c y4„„(k, cj), is calculated. Here, c is a normalization constant, \n) are the Kohn-Sham 
eigenstates, and 

^ , 2 -Im (nk|S(a;)|nk) 

A„„(k,w) = ^ — — 9 • (4) 

^ ' 7r[a;-£k-Re(nk|S(a;)|nk)]' + [Im(nk|S(u;)|nk)]' ^' 

To simulate the measured angle-resolved photoemission spectra from the calculated spec- 
tral functions, we multiplied the spectral function by the Fermi-Dirac distribution /(w) 
[Eq. (1) of the manuscript] with T = 25 K at which the experiments were performed^^ii^. 
Then, to take into account the experimental resolutions in energy and momentum, we have 
convoluted the intensity maps with a two-dimensional Lorentzian mask with Ak=Q.Ql 
and Aw = 25 meV, corresponding to the experimental resolutionii^i^. (This convolution 
results in finite linewidths even for zero binding-energy states.) In simulating the pho- 
toemission spectra along the FKM direction (Fig. 4a of the manuscript), we have used 
only one branch of the two linear bands in order to simulate the matrix element effects in 
/o(k, e,^, -Ej,) [Eq. (1) of the manuscript] (cf. Fig. 2 of Ref.—). In calculating the width of 
momentum distribution curves (the linewidths in Fig. 4c), we have arbitrarily subtracted 
off a constant from the simulated widths so that the width vanishes at zero binding energy, 
as also done in the analysis of experimental data in Ref.—. 

3. Substrate screening 

To include the effects of the dielectric screening associated with the silicon carbide (SiC) sub- 
strate, we have, as done in previous studiesii'^^'i^ (Supplementary Discussion 2 and Fig. |6]), 
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renormalized the bare Coulomb interaction by an effective background dielectric constant 
eb = (1 + esic)/2, where esic=6.6 is the optical dielectric constant of silicon carbide^^; one 
takes the average of the vacuum dielectric constant and the substrate dielectric constant 
because graphene is sandwiched in between the two media^^. Along the lines of Ref.--, we 
do not take into account atomistic interactions between graphene and the reconstructed 
surface of the silicon carbide substrate^^""— . 



Supplementary information 

1. Velocity renormalization 

Our simulated spectral functions allow us to study the velocity of Dirac fermions in the linear 
regime (away from the Dirac point) as a function of dielectric screening of the substrate and 
doping (Fig. E]). The band velocity of our model epitaxial graphene [eb = (1 + esic)/2 = 3.8] 
is found to be smaller than that of suspended graphene (eb = 1) by as much as ~10% 
(Fig. |5]). In addition, the velocity decreases as doping increases (Fig. |5]), in agreement 
with the previous calculation"^. Both trends are easily explained by observing that the 
polarizability of the substrate and the additional electrons in the graphene layer both lead 
to weaker e-e interactions. 

2. Substrate dielectric function 

In this section, we show that the dielectric function of silicon carbide (SiC) £:^'*"(q, is well 
represented by the value at q = and u = {i.e. , the optical dielectric constant) as far 
as our calculation is concerned. We also estimate that an error in the electron self energy 
arising from this simplification is less than 10%. 

In calculating the imaginary part of the self-energy, the relevant energy scale of the 
dielectric function of SiC is the quasiparticle energy measured from the Fermi surface, since 
the hfetime is determined by the real decay processes to lower energy states [see, e. g. , Eq. (5) 



of Ref. |35||. The energy argument of the inverse dielectric function used in the calculation 
of the imaginary part of the self energy at energy u varies between u itself and the Fermi 
energy due to the two Heaviside functions. As long as the dielectric function is reasonably 
constant over this range, the results are valid. Since we are only interested in the value of 
the linewidth for states from the Fermi level down to about 2.5 eV below it (as measured 
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by experiment and presented in Figs. 2 and 4 of the main manuscript), we have checked 
the vahdity of our approximation on substrate screening within this energy range and the 
corresponding wavector range as discussed later. The real part of the self energy is affected 
by the dielectric screening involving large wavevectors; however, the contribution coming 
from larger wavevector scatterings is smaller because the Coulomb interaction decreases 
with q. The error in the real part of the self energy arising from the inaccuracy in the 
dielectric function for a larger wavevector would be smaller than that involving a smaller 
wavevector, which is estimated below. 

In order to give a quantitative estimate of the possible error arising from the frequency and 
wavevector dependence of the dielectric function, we performed first-principles calculations 
of the dielectric function of 3C-SiC showing dielectric responses very similar to 6H-SiO^, 
which is the substrate used in experiments^iii^. (The unit cell of 3C-SiC is much smaller 
than that of 6H-SiC though.) 

Supplementary Fig. [6] shows that using the optical dielectric constant (w = 0, q = 0) is 
a good approximation for the energy and wavevector regime considered in our work: the 
maximum variation in the dielectric function is 15%. (We have also checked that the inverse 
dielectric function shows similar behaviors.) Moreover, since (i) the self energy is an average 
of many contributions, and, (ii) the finite frequency and finite wavevector effects increase 
and decrease the value of the dielectric function, respectively, the combined error coming 
from our approximation will be less than 10%. 

Supplementary Fig. [6] also shows (difference between solid and dashed lines) that in fact 
the maximum anisotropy in the dielectric function of SiC is ~3%. (We have also checked 
that the inverse dielectric function shows similar behaviors.) Since the self energy is an 
average of all contributions, we expect that the anisotropy in the calculated electron self 
energy arising from that in the dielectric function of SiC will be even smaller than that. 

3. Comments on other calculations 

First-principles calculations of the real part of the self-energy in graphene arising from e-e 
interactions within the G approximation have been reported previously^"^. The authors of 



Ref. 



J calculated the frequency dependence of the dielectric matrices within the generalized 



plasmon-pole model. The authors of Ref. [s] calculated the full frequency dependence of the 
dielectric matrices using the random phase approximation as we did for the present work. 
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For consistency we compare our calculations to the latter study. 

Our calculated velocity (1.23x10^ m/s) in suspended graphene is ~9% larger than the 
one reported in Ref. Isl (1.12x10^ m/s). Moreover, unlike the finding of Ref. Isi, we do not 
observe a kink in the quasiparticle band structure of suspended pristine graphene at an 
energy uo ~ —0.15 eV when e-ph interactions are not included. The fact that we observe a 
gradual increase in the band velocity when approaching the Dirac point energy, rather than 
a kink at a finite energy value is in line with results of model analytical calculations^i^I. 
To clarify this difference, we observe that, in the case of graphene, electronic states with 
wavevectors on nearest neighboring grid points of a discretized N x N mesh of the full 
Brillouin zone have an energy difference AE (eV) ~ 20/N. The Brillouin zone sampling 
adopted in Ref. 1^ (10 x 10 points) corresponds to electronic eigenstates with minimal energy 
separation of approximately 2 eV. We find that this energy resolution is not sufficient to 
achieve convergence in the dielectric matrices needed for the GW self energy (we used instead 
a 72 X 72 grid, corresponding to energy separations of about 0.25 eV). Therefore we suggest 



that t 
Ref. 



le difference between the band dispersions and velocities calculated here and those of 
may arise from the insufficient Brillouin zone sampling adopted in that work. 
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FIG. 1: Diagrams included in the calculated electron self-energy S. a, Electron self- 
energy S*^^^ = iGoWo arising from e-e interactions within the GqWq approximation. Go is the 
Green's function for bare electrons and Wq is the screened Coulomb interaction, b, Self-energy 
j^e-ph _ ig'^Q^D arising from e-ph interactions within the Migdal approximation. Here, g is the 
e-ph interaction matrix element and D is the dressed phonon propagator. 

FIG. 2: Simulated ARPES spectra, energy distribution curves (EDCs) and quasipar- 
ticle band structures of suspended graphene including e-e and e-ph interactions, a, 

Simulated ARPES spectrum of pristine graphene at T = 25 K taken along the Brillouin zone 
segment indicated in the inset of c. b, EDCs extracted from a. The central red curve corresponds 
to A; = (the K point), c, Quasiparticle band structure (solid red curve) obtained by connecting 
the peak positions of EDCs in b. d to f, and g to i, Same quantities as in a to c for n-doped 
graphene with charge densities of 4.5 x 10^^ cm~^ and 1.2 x 10^^ cm~^, respectively. The dashed 
blue lines in c, f, and i indicate the asymptotes of the linear bands far from the Dirac point. The 
energy difference between the upper and the lower asymptotes close to the Dirac point is indicated 
by Akink- 

FIG. 4: Momentum distribution curves (MDCs) of graphene and associated linewidths. 

a, Simulated ARPES spectrum of suspended n-doped graphene, for a doping level corresponding 
to a charge density of 4.5 x 10^^ cm~^, taken along the Brillouin zone segment indicated in the 
inset of c. b, MDCs obtained from a. c, Width of the MDCs obtained from b for suspended 
graphene (red curve) and that for graphene with a model dielectric screening corresponding to the 
SiC substrate (blue curve). The measured widths of the MDCsii are shown for comparison (black 
curve) . 

FIG. 3: Mismatch between the upper and lower bands of the Dirac cone. Calculated 
energy difference Akink between the asymptotic lines close to the Dirac point of the upper and the 
lower linear bands vs. doping for suspended graphene (red squares) and for graphene with a model 
dielectric screening (blue circles) corresponding to the SiC substrate (see Methods). The lines are 
a guide to the eye. 
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FIG. 5: Supplementary Fig. 1: Quasiparticle velocity in graphene. The slopes (quasi- 
particle velocity) of the linear bands far from the Dirac point (indicated by the blue dashed lines 
in Figs. Ic, If, and li of the manuscript) vs. doping. Red squares and blue circles are calculated 
quantities for suspended graphene and for graphene with a model silicon carbide (SiC) dielectric 
screening, respectively. Black triangles are DFT results within the LDA. The lines are a guide to 
the eye. 



FIG. 6: Supplementary Fig. 2: Calculated dielectric functions of SiC versus energy. 

Macroscopic dielectric function of 3C-SiC eQ'Q'(q, w) versus energy oj. Quantities for different 
wavevectors are shown in different colors. Solid lines and dashed lines show quantities for the 
wavevector q in two representative directions in wavevector space: F-X and F-L, respectively. 
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